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1. Introduction 

Cellular automaton (CA) can be used as discrete mathematical models for a variety 
of phenomenafl]. One remarkable feature for CA being used as a model lies in the 
fact that the setup is mathematically simple. Since both dependent and independent 
variables are discrete and a set of values which the dependent variable takes is finite, 
a rule table is more suitable for representing logical expression in a CA model. Thus, 
CA is free from mathematical conditions such as continuity and differentiability. 

Due to the complete discreteness of CA, mathematical tools for analyzing 
its solutions of CA are very different from those for analyzing solutions of other 
mathematical objects such as differential equation, finite-difference equation and 
coupled map lattice. Pattern analysis, graph theory and statistical mathematics are 
often used to analyze CA [2]. However, in seeking for a general solution starting from 
arbitrary initial data of a CA, this difference often becomes an obstacle. Compact 
expression of solution using special functions is rarely available for CA due to its 
discreteness. 

One of the bridges between discrete and continuous mathematics to solve this 
problem was found in the area of integrable system. By applying the following simple 
formula for the transformation of variables, 

lim £ log(e^/^ + e^/^) = max(^, B), (1) 

e-!-+0 

we can obtain CA from a finite-difference equation with continuous variables. This 
procedure is called 'ultradiscretization'Oll]. Solution of the finite-difference equation 
can be ultradiscretized by the same procedure. For example, we obtain a binary CA 
called 'box and ball system' from the discrete Lotka-Volterra equation through the 
above transformation[31 [51 [SI [7]. Moreover, the elementary CA of rule number 184 
(ECA184) is obtained from the discrete Burgers equation^. ECA184 is considered 
in the subsequent section. Note that we can obtain the differential equation and its 
corresponding solution by taking a continuum limit for a finite-difference equation and 
its solution. The continuous analogue to the discrete Lotka-Volterra equation is the 
Korteweg-de Vries equation and that to the discrete Burgers equation is the Burgers 
equation. Therefore three different discrete levels of equations are directly related 
each other through two types of limiting procedure. 

Replacing basic operations is another way to obtain an ultradiscrete equation from 
a finite-difference one. Equation ([T]) shows that the addition of the finite-difference 
equation corresponds to the max operation of the ultradiscrete one. The multiplication 
and division correspond to + and — respectively since we have 

( ^}^o) ^ log(e'^^' X e^/^) ^A + B, ( hmj e log(e'^/7e^/^) =A-B. 

Though the subtraction is not well-defined for the ultradiscrete equation, we can 
automatically obtain the ultradiscrete equation and its solution through a procedure 
of replacing x and / by max, -I- and — respectively if the subtraction is not included 
explicitly in the finite-difference equation and its solution. 

The algebra defined by operations of max, -I- and — is called 'max-plus algebra' [9]. 
Binary operations such as AND, OR and NOT can be easily converted to max-plus 
operations. Since CA can be defined by binary operations, a rule table of CA is 
usually transformed into a max-plus equation. In contrast to a binary solution to the 
original CA, the solution to a max-plus equation is real- valued. In this sense, both the 
equation and the solution are extended to the real value range when we use max-plus 



Max-plus analysis on some binary particle systems 



3 



algebra. We should point out that it is important but difBcult to express a rule table 
of CA by an appropriate max-plus equation in order to grasp the behavior of solution 
analytically. 

In the present paper, we propose max-plus expression for rule table, general 
solution and asymptotic behavior of a class of evolutional binary CA. The sum of 
values of dependent variables at all spatial sites is conserved through time evolution 
in this class. ECA184 is included in this class. Such CA can be used as an important 
digital model for transportation or flow phenomena. In regard to the asymptotic 
behavior of solution for this type of CA, a phase transition can often be observed 
[lOl [TT| . In order to clarify the asymptotic behavior mathematically, finding a general 
solution for this type of CA is necessary. It is found that the max-plus expression is 
more suitable in finding such a general solution. 

The present paper is organized as follows. In Section[21 we define the class of CA as 
mentioned above and show examples of solution and observed asymptotic properties. 
In Section [21 by giving the max-plus expressions of evolution rules, we present general 
solutions of the CA and, further, analyze the asymptotic behavior based on general 
solutions. In Section|4j concluding remarks are given. In |Appendix A 



some convenient 



formulas about max and min operations are listed. In Appendix B[ we prove the 
solution to a specific CA. In [Appendix C[ we prove the asymptotic behavior of the 
solution. 



2. Classification of binary particle CA 

Let us consider the following evolution equation, 

uJ+i=/KVn •■-"J+rJ, (-oo<j-<^, 0<n) (2) 

where n and j are integers, representing the time step and the space site number, 
respectively, ri and r2 are also integer constants. Assume that the value of u takes 
either or 1, thus the above equation ([2]) gives an evolution rule of binary CA. The 
binary function / has R {— r2 — ri + 1) arguments and we call this CA 'i?-neighbors 
CA'. 



2.1. Identification of CA 

Using a transformation of coordinates defined by j ^ j — cn, the above equation ([2]) 
can be rewritten as 



u 



Thus, the solutions to these two equations ((2|) and Q can be transformed from one 
to the other. Therefore, we can identify equation ^ for any c with ([2]). 
Applying the transformation j — >■ —j to ([2|) , we obtain 

«r'=/K-.„«--.,+i, J, (4) 

and can identify ((4|) with ([2]) by the same reasoning. Moreover, if we 'flip' the 
dependent variable m, that is, if we apply the transformation u — !■ 1 — u, it follows that 



u 



n+l 
j 



Although the function g is generally different from /, we can still identify ([5]) with 
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Moreover, there are cases where f oi ^ does not depend on ah of its arguments, 
for example, in the foUowing equation 

u:;+'^f{u],u]^„u]^,), f{a,b,c)^g{a,b), 

the number of formal arguments of / is three (3-neighbors) but / does not depend on 
the last argument. Therefore, this equation is reduced to m"^^ — which 
is of 2-neighbors CA. Let us consider another equation, 

in which, the values of u's at even (odd) sites at time n + 1 depend only on those 
at even (odd) sites at time n. Thus, this equation is reduced to u"'''^ = g{u^ ,u^_^_i), 
which becomes 2-neighbors CA by the separation of even and odd sites. In the present 
paper, we exclude the equations which are reductive to those of less than i?-neighbors 
from a class of i?-neighbors. 

2.2. Particle CA 

Let us consider the following condition for any n, 

oo oo 
j — ~oo J — — oo 

which means the sum of u, that is, the number of I's at all space sites is conserved 
for arbitrary n. Let u" denote the number of particle at j-th site and n-th time 
step, and each '1' in the solution represent a particle. Therefore, the particles move 
among the sites according to the evolution rule defined by ^ without creation or 
annihilation. We call the CA satisfying the condition (jS]) 'particle CA' (PCA) in 
this sense. Classification and statistical analysis of PCA was reported by Fuks and 
BoccarafTUlfT^. 

If a CA defined by ([2]) is PCA under appropriate boundary conditions, it can 
always be written in a conservation form, 

where (7(01, 02, . . . , a_R-i) is a 'flux' for the 'density' ttjl3. ■ The condition 
g(0, 0, . . . , 0) = is necessary for a particle CA because a system with zero particle 
will produce zero flux. This condition can be used to determine an indefinite constant 
in the expression of q. 

In what follows, we list up all PCA of 1-neighbor to 4-neighbors identified in 
Subsection 12.11 Label 'n' or 'n-m' is assigned to each PCA where n denotes the 
number of neighbors, and m denotes the subclass of the same neighbors CA if it 
exists. Rules of PCA are defined in the form of ([7]) together with a rule table of q. 
The upper row of the table shows all possible combinations of binary arguments and 
the lower row gives the values of q. 

• PCAl: u^'+i g = 

• PCA2: none 

• PCA3: = + q{uj_^,u^) - g(Uj",u^+i) 



a b 


11 


10 


01 


00 


q{a,b) 





1 
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• PCA4: u 



4-1 



4-2 



4-3 



4-4 







111 


110 


101 


100 


oil 


010 


001 


000 


g(a,6,c) 





1 





1 





1 










abc 


111 


110 


101 


100 


oil 


010 


001 


000 


q{a,b,c) 





1 














-1 







abc 


111 


110 


101 


100 


oil 


010 


001 


000 


q{a,b,c) 





1 






















abc 


111 


110 


101 


100 


oil 


010 


001 


000 


q{a,b,c) 

















1 









For simplicity, let n = represent the initial time step of the evolution. For PCAl, 



we have always a trivial static solution 



PCA3 is equivalent to ECA184 [T], 



which is related to a simple traffic congestion model [8]. In what follows, we will 
discuss the dynamics of particles and the asymptotic behavior of solutions for PCA3 
and 4. Although the max-plus analysis of PC A3 (ECA184) was already given in [8], 
we here give a briefly review as a typical scenario of max-plus analysis. 

From the rule table of q, we can determine the dynamics of a particle system. 
In what follows, we show the motion representations of PCA3 and PCA4. Note that 
'1' denotes a site occupied by a particle and '0' denotes an empty site. The motion 
representation for each PCA gives the configurations in which only neighborhoods 
resulting in a particle motion assuming that by default, in all other cases, particles do 
not move 11 . The arrow indicates the initial and final positions for a moving particle 
in one time step. 



PCA3: 

PCA4-1: 

PCA4-3: 



10 
101 
110 



100 



PCA4-2: 
PCA4-4: 



110 
/> 
010 



001 



2.3. Asymptotic behavior of PCA 

In this subsection, we discuss the asymptotic behavior of solution from arbitrary initial 
data as n — oo . The behavior is mainly dependent on the density of particles, or the 
average number of particles per site. To observe this dependence, a periodic boundary 
condition for space sites is assumed. Let K be the number of sites in one period, the 
density p is defined by 

Since the number of particles is conserved, p is also a constant irrespective of time n. 

For most of PCA, the asymptotic behavior changes at certain critical values of 
p. The relation between the density and the average fiux per site reflects this critical 
phenomenon well. The average flux per site is deflned by 
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For large enough n, the evohition of a PCA often approaches a steady state, in which 
becomes a constant. In the case of convergence, the constant is defined by 

Q = hm 9". 



The vahie of Q is independent of the initial distribution of particles and depends 
solely on the density p for PCA3, 4-1, 4-2 and 4-3, whereas for PCA4-4, the value of Q 
depends on both the density p and the initial condition. The graph of Q vs. p is called 
'fundamental diagram' in the area of traffic flow analysis. We use this terminology 
throughout the present paper and abbreviate it as FD. The FD's for PCA3 and PCA4, 
as well as examples of evolutions are shown in Figures [T][5] respectively. 

The relations between the solutions of steady states at large enough n and the 
critical values of p are listed below. 



PCA3: u]- 



(0 < p < 1/2) 
(1/2<P< 1) 



PCA4-1: 



n+l 



(0 < P < 1/3) 

u]+, (1/3 < p < 1) 

<+i (0 < p < 1/3) 

PCA4-2: u]+^ = \ u]_^ (1/3 < p < 2/3) 

^J'+i (2/3 < p < 1) 

u] (0 < p < 1/2) 

PCA4-3: u'J+^ ^ I (1/2 <P< 2/3) 



PCA4-4: 



^^+1 (2/3 < p < 1) 

u] or u]_i (0 < p < 1/2) 
< (1/2 < P < 1) 



3. Max-plus analysis of PCA 

3.1. PC A3 0/ 

The flux function of PCA3 was found to be q{a, b) — min(a, 1 — 6) if a and b are binary. 
Thus, the evolution rule can be rewritten in the following max-plus form, 

uj+i = u"^ + min(u7_i, 1 - uj) - mixviu"^ , 1 - m^^i)- (») 

Using a transformation of variable from m to /, 

^-f?-fU + \^ (9) 
we obtain the evolution rule for / 

/;+i^max(/JLi,/;Vi)- (10) 
Note that, for simplicity, an arbitrary constant occurring in the process of deriving 
pllj) is set to 0. 

The initial value problem of (|T0| can be easily solved as 

/" — max(/j'_„, /j_„_|_2, . . . , /j'_|_„_2i fj+n)- 
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The inverse relation of © is 



fe=0 

with an arbitrary constant c, from which we have 

j+K 



-K - fj - X! 

k=j+l 



Therefore, it is obvious that /JY^r < /" for P < 1/2 and > /" for p> 1/2. 

Using above inequalities, we have 

max(/0_„, /j'^.^^^, ■ • • , f^-n+2K) (0 < P < 1/2) 

max(/jV„_2A:: f°+n-2K+2, fj+n) (1/2 < P < 1) 

for large enough n. Furthermore, we conclude that 

/;+^ = wj+'^u^^, (o<p<i/2) 

/r' =/"+!' "r'=^"+i (i/2<p<i) 

for large enough n. These analytical relations explain the steady solution of PCA3 
shown in Figure [1] 

Finally, substituting the relation u^^^ = under the case of0<p<l/2 into 
the evolution equation ([S]) , we arrive at 

- min(-u"_i, 1 - u") = u" - min(w", 1 - 

for any j. In other words, u"_i — min(u"_2i 1 ^ '^j) ^ constant irrespective of j. 
Since < p < 1/2, there exists at least one site j satisfying = (0,0) or 

(1,0), which resulting in a zero value for the above constant. Thus, we have 

in the case of < p < 1/2. Similarly we obtain Q = 1 ~ p in the case of 1/2 < p < 1. 
This relation gives FD of PCA3 as shown in Figure [1] 



3.2. PCA4-1 

This PCA is known as the Fukui-Ishibashi model used for traffic congestion 
analysis 1 14) and its Euler and Lagrange representation are studied using max-plus 
expression [15]. Asymptotic behavior for the solution of PCA4-1 is similar to that of 
PCA3 except a few slight differences: i) critical value of p is 1/3, ii) for the smaller 
region of p (0 < p < 1/3), all particles move with a speed of 2 at large enough n. The 
evolution equation of PCA4-1 can be described in the following max-plus form, 

= u" + min(u"_2 + 1 - u") - min(u"_i + u", 1 - u"+i). 

Introducing a transformation 

^-f?~fU + \^ (11) 

we obtain 



= max(/;Y2,/;Vi) 
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by choosing an appropriate value for the arbitrary constant. 

Consequently, the general solution for the initial value problem turns out to be 

■fj = ™9'^(/j-2nJ /f-2n+3' ' • ■ ' /j+n-3' fj+n)- 

On the other hand, the inverse relation of (|Tl]) gives 

fc=0 

which leads to 

f7+K-f7 = K{p-l 

By similar calculations for PCA3, at large enough n, we can derive steady state 
solutions: 

cn+l _ fn „,"+! 



/r =/"-2, =^-2 (o<p<i/3) 

/;+i = /jVi, u«+i = ^7+i (i/3<p<i) 

and further the following relations 



min(u"_i + u", 1 - = 



u]_i+u] (0<p<l/3) 
(l/3<p<l) 



From above relations, it can be easily shown that 

p^hP iO<P< 1/3) 

\l-p (l/3<p< 1) ■ 



(12) 



In summary, the above analytical results agree with the solutions and FD shown in 
Figure [H 

3.3. PCA4-2 

It can be shown that the evolution equation of PCA4-2 is 



u 



n+l 



(13) 



I q{a, b, c) = min(max(— c, a + & — 1), 1 — c) 
Through a transformation, 

u]^f-- /;u + 1, (14) 

it follows, from (IT51) . that the evolution equation to / is 

= max(min(/j^2,/;Vi + I) Jj+i (15) 
The general solution to equation (fT5|) turns out to be 

/" = max fminf/y_n,, I ni, — — , min /? c.,- + -1 V (16) 
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By expanding RHS of (IT5|) . specific examples of for n = 1, 2, 3 are given as 

= max(min(/^_2, /jVi + ^) . /jVi - ^) - 
/2 = max(min(/J'_4, ./^^i + ^ /jVa + l) , min(/0_i - i /jVs + l) , - l) , 
ff = max(min(/0_e, /J'^g + 1 /« + 1, f^^, + ^) , mm[f^_, _ 1 /« + 1, /j^g + ^) , 

min(/]'-l,/jV3 + ^),/jV3-^). 

The proof for the general solution (|16p satisfying ([T5|) is given in | Appendix B| 

Moreover, at large enough n, the asymptotic behavior of the solution, as well as 
the value of Q, can be shown to be: 

fo<.-<^-.(^^"-- + ^) (0<P<l/3) 



(17) 



L,f#K^-A^^"-- - (2/3<P<l) 

/;+' = /;Vi + ^> «r'-"5Vi (o<p<i/3) 
/;+^ = /jVi-i ""+' = ""+1 (2/3<P<i) 



= /"-2, = <-2 (1/3 < P < 2/3) , (18) 



2' 

r -"7 (0 < p < 1/3) 

= < ""-2 + <-i - 1 (1/3 <P< 2/3) , (19) 

I 1 - (2/3 < p < 1) 

r (0 < p < 1/3) 

Q = < 2p- 1 (1/3 <p< 2/3) . (20) 

[l-p (2/3<p<l) 

The proof of the above results is given in [Appendix C| In summary, the analytical 
results reveal the steady state solutions and FD as shown in Figure [31 

34. PCA4-3 

The evolution equation of PCA4-3 can be expressed by 

{n+l n , / n n n\ / n n n \ 

q{a, b, c) = min(max(0, a + 6 — 1), 1 — c) 
By using a transformation, 



rn r7l I 

Jj Jj-1 + X, 
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we arrive at an evolution equation for /, 

= max(min(/j^2, /;), /jVi - ^) • (21) 
From which, the general solution turns out to be 

/" = max fminff? n„ , n., - min /?_2.,4-9,D 

0<k<n\ ^"+'^'' 2 3k/2<i<n^ 2n+2iJJ 

by means of a similar proof as PCA4-2. Therefore, it is omitted here. 
The specific expressions of /" for n = 1, 2, 3 are as follows: 



/j=max(min(/^_2,/]'),/]Vi 



ff = max(min(/^_4, /j'^^, /]>), min(/^_i - 1, /]>) , /« 
/3 = max(min(/^_e, /j'^^, /«), min(/0_3 - 1, ,f°_2, /°) /jV 



3^ 

At large enough n, we can derive the asymptotic behavior of solution and, further, 
the value of Q, which both depend on the density p as follows 



(0<P<l/2) 

^max ,f minf/°_2„+3fe 



o<k<K-i\ yj-^n+.i^ 2' 

1<?<K-1 ■'j-2»+2L3fc/2j+2t 



(1/2 <p< 2/3) 



n 



- k 



Lo,f#K^-A^^"---V'j (2/3<p<l) 

= - (0 < p < 1/2) 

/r+' = = ""-2 (1/2 < p < 2/3) _ 

2' 



iT^'-/;Vi-^> ^r'-"^! (2/3<P<i) 



(0 < p < 1/2) 
q{u]^„u]_„u]) ^ { u]_, + u]_, - 1 (1/2 <p< 2/3) , 

1 - (2/3 < p < 1) 
'0 (0<p<l/3) 

Q = I 2p-l (1/3 <p< 2/3) . 
_l-p (2/3<p<l) 

Again, the above results reveal the solutions and FD as shown in Figure HI Since the 
proof is tedious and similar to that of PCA4-2, it is omitted here. 



3.5. PCA4-4 

We are unable to find a general solution to this type of PCA in the max-plus form. 
However, it is worth pointing out the following observations from the binary rule table 
in the case of0<p<l/2: if there exists a site j satisfying u^_^ = u° = 1 in the 
initial data, the asymptotic solution will become static, that is, u^^^ = u" at large 
enough n. If there is no two consecutive I's in the initial data, the solution satisfies 
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u^_i — for any j and n. It then foUows that Q = in the former case and Q — p 
in the latter case. Note that two different types of initial data are recognized, which 
correspond to two branches in FD as shown in Figure [S] in the case oi p < 1/2. 



4. Concluding remarks 

We propose the max-plus expressions to equations and solutions of PCA4-1, 4-2 and 
4-3. We analyze the asymptotic behaviors of solutions and prove the fundamental 
diagrams mathematically. The max and min formulas are effectively utilized in the 
proofs. In view of the role of basic operations appearing in max-plus algebra, the way 
of proofs is similar to that for the finite-difference or differential equation. In this 
sense, the mathematical technique used is analytic. It is important to introduce such 
an analytic technique into the digital systems since it provides us with an alternative 
approach to study them, and sometimes can lead to rigorous results. 

We restrict ourselves to the condition that the value of u is binary in the proof of 
solutions for PCA4-2 and 4-3. Hence, the max-plus expressions for general solutions 
proposed only hold under this condition. On the other hand, solutions for PCA3 
and 4-1 are free from this condition. Their solutions and evaluation of asymptotic 
behaviors from more general data (real valued data) are valid. It is a future problem 
to check the validity of the present expressions for PCA4-2 and 4-3 from more general 
data. 

Finally, we have not yet obtained the max-plus expression for solution to PCA4- 
4 and PCA of more than four neighbors. The general relation between PCA and 
its corresponding continuous equation is not clear by using the ultradiscretization 
method. It is worthy to solve these problems in the future. 



Appendix A. Formulas about m£LX and min operations 

There are useful formulas to derive or prove relations concerning max and min 
operations. Below we list up such formulas relevant to this paper. 
The definition of max and min operations are as follows. 

(a if a > &, . / , X f " if a < 5, 

• max(o, o) = < , mm(a, oj = < . 

I b otherwise I b otherwise 

From this definition, max and min are related each other. 

• max(a, 6) = — min(— a, —6), min(a, b) = — max(— a, —6). 

Most basic formulas satisfy commutative, associative and distributive laws as follows, 
where max or min plays a role of the addition and i.e., the usual sum of real 
numbers serves as the multiplication. 

• max(a, b) — max(6, a), min(a, b) = min(6, a). 

• max(max(a, 6), c) = max(a, max(6, c)), min(min(a, b),c) — min(a, min(6, c)). 

• max(a, b) + c = max(a + c, b + c) , min(a. b) + c ^ min(a + c, b + c) . 

From the commutative and associative laws, we can easily extend max and min 
operations from binary to multinary. 

• max(ai, 02, . . •) = the maximum value among ai, 02, . . ., 
min(ai, 02, ■ ■ ■) — the minimum value among ai, 02, . . .. 
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We have the following formulas about these multinary operations. 

• max(ai, 02, ...) = — min(— ai, —02, . . .), 
inin(ai, 02, . • .) = — max(— ai, —02, . . .). 

• max(. . . ,ai, . . . ,aj, . . .) = max(. . . , , . . . , a^, . . .), 

min(. . . ,a.i, . . . ,aj, . . .) = min(. . . , Oj, . . . , a^, . . .). 

• max(niax(ai, 02, . . .), ^2, • • •) = niax(ai, 02, . . . , 61, 62, • • 
min(min(ai, 02, . . .), 61, &2, • ■ •) = min(ai, 02, . . . , 61, 62, . . .)• 

• max(ai, 02, . . .) + 6 = max(ai + 6, 02 + &, . . .)i 
min(ai, 02, . . .) + b = min(ai + b,a2 + b, . . .). 

• If there exist i such that a < bi, max(a, 61, 62, • ■ •) — niax(6i, 62, • • ■)■ 
If there exist i such that a > bi, min(a, 61, 62, ■ • •) — niin(6i, 62, • • •)■ 

• max(min(ai, 02, . . ^1,^2, ■ ■ ■) — min(max(ai, &i, 62, ■ • ■), inax(a2, &i, &2, ■■■),■■■) 
min(max(ai, 02, . . bi,b2, ■ ■ ■) = niax(min(ai, 61, 62, . . .),uim{a2, 61, 62, ■■■),■■■) 



The last set of formula regarding the exchange of max and min operations has been 
widely used in this paper. 

Appendix B. Proof that (fT6|) satisfies ([TSj) 



We prove that (|16p satisfies (|T5|) by the mathematical induction. For n = 1, it is true 
since (fT6|) is equivalent to (fTSj) of n = 0. Let us assume that (fT6|) satisfies (fT5|) for 



n — m (m > 1). Then is expressed by fj as 



Using we obtain 



= max min 



maxfminf /°_2m+3/£-2J /j-2m+3fc+l + 9) ■ fj-27n+3k+l 



2) 2' 

1\ „n 1\\ i 



fc+]Xi<m(™'^^(™^^(-^^-2m+3i-2>/i-2m+3»+l + ^ J^-2,n+S.i+l ^2)) ^2 
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Using the max-min formula, we have 

■> 0<fe<m V 



min^max^/°_2TO+3fe_2, /j_2m+3/s+i ~ 2)' 

inax^/°_2„_|_3fc_|_i + 2' /j'-2m+3fc+l ~~ 2 
fe+]xi<m(™'^(™^^('^^°-2'"+3»-2'4'-2m+3i+l " ^) . 

max^/°_2„+3j+i + 2'/j-2m+3j+i "2))) 



= max mm 

0<fe<m V 



max^/j_2jjj_|_3fc_2 2'/7-2m+3fe+l 2 ) ' /7'-2m+3fc+l 2 ' 

fe+l<i<m(™^'''('^'''~^"*"'"^'~^ 2 ' •^•j°-2m+3i+l + 2 ) ' •^i-2m+3j+l + 2 ) 



Expanding the terms, we have 
= max fminf 

0<fe<m\ V 

Jj-2m+3k-2 ~ 2 ' ■'i-2m+3fe+l ^ J > /j 

inax^/j_2TO_l_3fc_l_l H ^ J /j-2m+3/c+4 2 / ' -'^-2^1+3/0+4 



i-2m+3fc+l ^5 • 



max^/j_2TO+3fe+4 -I- 2 ; /j-2m+3/c+7 + 2 y ' j-2rn+3fe+7 + 



mdX^_/^_l_^_2 -|- 2 ; J j+m+1 2 J'Jj+m+1^ ^ 



2 ' 



Since the underlined terms are neghgible, we obtain 

A; ,n fc + 1 ^ 



/f^' = o^fc^„(min(max(/°_2„+3fe_2 - 2./°-2m+3fc+i 



2 



m+iXN 

■•'Jj+m+lT^ 2 7/ 



Jj-2m+3k+l 2 ' Jj-2m+3k+4 2 ' •' j-2m+3ft+7 

TO + 1 ' 
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Using the max-min formula again, we have 



2 ' 

J j~2m+3k+4 ' 2 ' J j-'im+Zk+l 'i 2 ' ' • ' ' Jj+m+1 ' 2 



-2m+3/c+l 



/C+ 1 



J j-2m+3k+4 + 2 ' •'j-2m+3fc+7 ' 2 ' ' " ' ' J j+m+1 ' 2 

Expanding the terms, we have 

fj ^ = max^ min^/j'_2m-27 fj-2m+l + 2 ' fj-2m+4 + 1) • • • I /j+m+l ^ ^ ^ i 

minffO -if" +1 f" +- f" +!!i±I'| 

minff" -if" f" +- f" +!!i±l'| 

™'^(-^i-2m+4 ^ 1' /i-2m+7 + 2' /i-2m+10 + 2, ■ ■ • , /j+m+l H ^ 
min^/"_2m+4 ~ 1' fj-2m+7 ^ 2 ' ■fj-2m+10 + 2, ■ • ■ , fj+rn+l ^ ^ 

™-i^(/i-2m+7 ~ 2' /f-2m+10 + 2, /]'-2m+13 + 2' ' ' ' ' -^i+m+l ^ 2 
. . . , 

. / _ m-1 _ m - 2 m+l 

™.UH^/j_|_m_5 2 1 /j+m-2 2 ' •fj+m + 1 + 2 

TO TO + 1 

mm ' ' 



\^''i+m-2 2 ' 2 / 

min|/j^jjj_2 2 ' -fi+m+i 2 ' -fi+ni+i ' 2 



■ / .0 ™+ ^ .0 , m + i\\ 

miiH^/j+„j+i ^ 7 /j+m+i "I 2 J ) ' 

Since the underhned terms are neghgible, we obtain 

fj"^^ — max^ min(^/j_2m-2' /j°-2m+l + 2 ' /]'-2m+4 + 1; ■ • • I /j+rn+l ^ ^ ) J 



min^/"_2m_|_l 2 ' /i-2?n+4 + 1; fj-2ni+7 + 2 ' ' ' ' ' ^ 2 



lin^/"_2m+4 ~ fj-2m+7 + 2 ' /i-2m+10 + 2 



3 .n , „ .0 , m+l 



. /,0 .0 TO + 1 

■'J+m-2 ' /j+m+1 T ^ 

TO + 1 ' 



_TO^\\ 
./ j+m+1 2 // 

( min(/"_2„+3fc-2 - min /"_2 

L\ \ 2 fe+l<i<m+l 



0<fc<m+l V ■ yj-^'u-r^^-. 2 ' k+l<l<m+l ■'j-2m+3^-2 + 2. 

The above equation is nothing but (jl6p for n = to + 1. Thus the proof is complete by 
the mathematical induction. 
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Appendix C. Proof about relations from (|17p to (|20p 

The asymptotic behavior of solutions to PCA4-2 is different for three ranges of p, 
< p < 1/3, 1/3 < p < 2/3 and 2/3 < p < 1. Therefore we consider three cases 
separately. First, we expand the solution (ITBl) for rt » as follows to make evaluation 
easier; 

/j" = max^ niill(/j'-2ni /j'-2ri+3 + 2 ' fj-2n+6 + 1; ■ • ■ J fj+n + 2 

min^/j'_2„+3 ^ 2' ■^j"""2n+6 + 1' fj-2n+9 + 2' ' " ' ' ^ '2 

min^/j'_2„+6 ~ 1; /j-2n+9 + 2 ' /j°-2n+12 + 2, ■ ■ ■ , + ' j-Q 2) 

. . . , 

mm(^/j+„_3 —Jj+n + 2 



Jj+n 



Appendix C.l. Case of < p < 1 /3 
In this case, we have 



j+3pK 



for any j and 71 where p is a non-negative integer and K is the spatial period. Therefore 
we have 

/;>/;V3pK + ^- (C.2) 

Let us define by 

g" = min f/?-i H — 

■ f fO n-K + l ri-X + 2 g n\ 

— mm^/j_|_„_3;f^3 i ^ , Jj+n-3K+6 ^ ^ ' ■ • ■ ' Jj+n ^2)' 

For p > and < q < K, 

ro , n-pK-q - g 

/j+n-3(p/f+9) + 2 - J+n-39 + 2 

holds from (IC.2p . Therefore the first min term of (jC.l[) is equal to for n ^ 0. We 
can make a similar evaluation from the second to the (n — K + l)-th min terms of 
(fell) as 

inin^/j'_2ii+3 ^ 2 ' -^j-an+e + 1' fj-2n+9 + 2 ' ' ' ' ' -^5°+" ^ '2 ) ^ ™l^(/j-ri+3 ^ 2 ' ) ' 
(/i-2ji+6 ^ li fj-2n+9 + 2 ' •^j-2n+12 + 2, • • ■ , fj+n ^ '2 ) ^ niin^/°_2„+6 ^ .9j ) i 



mm 



innH^/j+„_3A' - J Jj+n-3K+3 ^ ^ ' J]+n- 



n- K ,n n-X + l.n n-X + 2 

-3K+6 "I ;^ ' 



Q n\ ■ f rO n ~ K 
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Moreover, the relation 

holds since w" = or 1. It implies that the difference between /" and fj±p is finite 
if p is finite. Therefore we can evaluate from the {n — K + 2)-th to the penultimate 
terms of (jC.l|) for n as 

• frQ n-K + 1 , n-K + 2 g n\ 

mni^/j+„-3_ff+3 ^ . Jj+n~3K+6 ^ ^ ' • ■ • ' Jj+n ~^ 2 J 

_ n-K+1 

— J j+n-3K+3 2 ' 

• frO n-K + 2 , n-K + 3 „ n\ 

mm yJj+n-3K+e ^ > ij+n-Sif+g ~1 ^ ' • • • ' /j+n 2' ^ 

n — 2 „n JT- — 1 „n ?J\ „n n — 2 



M^Jj+n-e 2 ' j+n-3 I 2 'Jj+n' 2) ^ 

mmi /j^„_3 - , + ) — /j+n-3 7^ ■ 



y]^n-^ 2 '■'j+'i ' 2, 

Thus (|C.1|) becomes 

= max(g^", min(4^_„+3 - i, s}') , . . . , min(/°+„_3^. - 5^") , 



w- JiT+l 

jj+Ti-s-ff+a 2 ' ■ • ■ ' •/j+n 2 



It coincides with (|17l) in the case of < p < 1/3. 
Since we can easily show 

5j -5j+i + 2' 

we have 

1 ,,1 
2 

Substituting w"^^ = ""j+i m^*^ evolution equation (jl3p . we obtain 

Therefore w" + g(u"_2, w") is a constant irrespective of j. Since there exists at 
least a site j such that (u"_2, w") = (0,0,0) or (1,0,0) for < p < 1/3, we 
obtain m" + g(-u"„2i "j-ii ) = ^'i^Y J- Thus we have 

Q = :^ E "7) = E - -p- 

The above results coincide with (fT8)) . (fT9|) and ((20l) in the case of < p < 1/3. 



Max-plus analysis on some binary particle systems 17 

Appendix G.2. Case ofl/3<p<2/3 
In this case, the inequahty 

fj^+3pK ^ ^ /j" ^ fj+3pK ^ ^> (C-3) 

holds for any j and n where p > 0. Under this condition, the first to the {K — l)-th 
min terms of (jC.ip for n ^ are evaluated as follows: 

min(/«_2„,4'_2„+3 + ^:---:/]Vn + f) 

= min|^/°_2„! fj-2n+3 + 2' • ■ • ' fj-2n+3{K-l) ^ ) ' 

min^/j'_2ii+3 ^ 2' '^j-an+e + 1: • ■ • j + '^^ 

= niin^/°_2n+3 ^ 2' •^l'-2n+6 + • ■ • i fj~2n+3K + ' 



™^'^l^''i-2n+3(if-2) 2 ' ''j-2n+3(if-l) 2 ' ' ' ' ' '2/ 



— niin|^/^_2„_|_3(^_2) 2 ' ''j-2n+3(A:-l) ' 2 ' 

2ii: - 3n 



.0 +f£l_J^^ 

• ■ • ' Jj~2n+3(2K-3) ^ 2 y ' 

Regarding the rest min terms of (|C.ip . they are expressed generally by 

r = mm /,_2«+3fc " TT' . ^in / _2' 



Z k+l<t<n Z/ 

where K — 1 < k < n. We can evaluate them similarly as above and obtain 

rk = min(/0_2„+3fc " |' ... /°-2n+3« + ^) ■ 

Then the inequality 

j-2n+3i , 

Q 1\ I + A: 



(/j-2n+3i + 2) (/j-2n+3fc 2) ~ 



=j-2n+3fc+l 



2 



> ^— — - H — = 2k-i>2k~ min(n, fc + JiT - 1) 

= max(2fc - n,k - K +1) >0 

holds since i < min(n, k + K — 1) and k > K — 1. Thus we obtain 

- f ^ 

^fe — /j-2ri+3fc — 2' 

Moreover, note that 
holds for p > by (IC31) . 



It can be easily shown that the relation /"^^ ~ 17-2 holds for /" in the above 
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Therefore we obtain the final results about the asymptotic evaluation of /J' for 
n ^ in the case ofl/3</3<2/3as 

/" = max^ min^/j'_2„, /]'-2n+3 + 2 ' ' • ■ ' ^°-2n+3(/f-i) ^ 2 
inin^/j'_2n+3 ^ 2^ f^-2n+& + 1: • ■ • 1 fj~2n+3K + 

• ffO fO , K-l 

f +^-^) 
■ ■ • ' J j-2n+3{2K-3) ' 2 / ' 

fj-2n+3{K-l) 2 ^fj-2n+3K ■ • • J fj-2n+3{2K -2) ^ ^ 1)) • 

It coincides with (|17l) in the case of 1/3 < p < 2/3. 
equation. Hence so does the relation u"^^ = uj_2- Substituting it into (fT3|) . we obtain 

n , n in n n\ n , n in n n \ 

"j-2 + "i-l - q{Uj-2,Uj_i,Uj ) = Uj_^ + - q{u^_^,Uj , Uj + i), 

and u"_2 + — (7(u"_2' "i-n ) ^ constant irrespective of j. Since there exists 
at least a site j such that (m"_2, u^-i,u'j ) = (1, 0, 0) or (1, 1, 0) for 1/3 < p < 2/3, we 
have u"_2 + — (7(u"_2, m") 1. Thus we derive 

Q=^f2 '?K-2, ^-1, "7) = ^ EK-2 + ""-1 - 1) = 2p - 1. 

The above resuhs coincide with ^ and (UHl) in the case of 1/3 < p < 2/3. 
Appendix C.3. Case oj 2 / 3 < p < 1 

Substituting the transformation u — 1 ~ u into (|13p . we obtain 

I q(a, 6, c) = — min(max(c — 1, 1 — a — fo), c) 
Using the max-min formulas, we have 

5(0, h, c) = max(min(l — c, a + & — 1), — c) = min(max(l — c, — c), max(a + 6 — 1, — c)) 
= min(max(— c, a + & — 1), 1 — c). 

Since q is the same as g, the evolution equation ([T3| is self-conjugate, that is, it is 
invariant under the above transformation. Therefore, the behaviors for solutions of 
< p < 1/3 and 2/3 < p < 1 are symmetric to each other and there exist one-to-one 
correspondence between the solutions. Since u"^^ — uj_^i holds for n ^ in the case 

of < p < 1/3, 1 — w"^^ = 1 — that is, m"^"^ = "j+i ^g^in holds in the case 

of 2/3 < p < 1. Substituting this relation to p^ . u" -I- g(u"_2, w") is again a 
constant. Since there exists at least a site j such that (u"_2, m") = (1, 1,0) or 
(1, 1, 1) for 2/3 < p < 1, we obtain u" +q(w"_2 7 ""j-ij ) = 1 foi' ^^y j. Thus we have 
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These results coincide with ([H]), (HI]) and 1^ in the case of 2/3 < p < 1. 

QS 
1 



Moreover, introducing the transformation /" — — /" with = 1 — w", we obtain 



from (HU, and 

= min(max(/j^„ /;Vi - ^) ■ /JVi + ^ 

= max(min(/;U,/;Vi + ^)'/;-i-£ 

from . These relations are in the same form of (ITil) and p3|) . Thus the asymptotic 
solution /" for 2/3</9<lis derived from that for < p < 1/3 as 

It coincides with (|17p in the case of 2/3 < p < 1. 

Note that the above symmetry of bit inversion holds for any n as well as for n ^ 
and it gives another expression of solution as follows: 

/" = — max fminf — /°_o„ , — — , min — . ■h, + 



min fmaxf /? n„ I ^i, + — , max /°_ 



k „r\ i 

2n+3z 7y 
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Figure 3. FD and solutions of PCA4-2 
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Figure 4. FD and solutions of PCA4-3 
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Figure 5. FD and solutions of PCA4-4 



